ABSTRACT
Cellulose synthase (CesA), one of the key enzymes in the biosynthesis of cellulose in plants, plays an important role in plant growth and plant resistance. In this study, a total of 21 AsCesA genes from Aquilaria sinensis were systematically identified and the physico-chemical characteristics were analyzed based on genome database and bioinformatical methods. The phylogenetic tree was constructed and the gene location on chromosome, cis-acting elements in the 2 000 basepairs upstream regulatory regions and conservative motifs were analyzed. The AsCesA proteins were mainly located on the plasma membrane. The number of amino acids of the proteins ranged from 390 to 1 261. The isoelectric point distributed from 5.67 to 8.86. All of the 21 AsCesA proteins possessed the transmembrane domains, the number of which was from 6 to 8. The genes were classified into 3 groups according to the phylogenetic relationship. Obvious differences were observed in motif composition in genes from different groups. However, motif2, motif6, motif7 and motif10 were observed in all of AsCesA proteins. Analysis of cis-acting elements indicated that AsCesA genes family has cis-acting elements related to plant hormones, abiotic stresses, and biological processes. Seven AsCesA genes with differential expression were selected according to the calli transcriptome data induced by NaCl at different times and their expression levels under different abiotic stresses were analyzed by quantitative real-time PCR. The results indicated that salt, low temperature, drought, and heavy metal stresses could affect the expression level of AsCesA genes, and the abundance of AsCesA1, AsCesA3 and AsCesA20 showed a significant change, implying their potential important roles to the abiotic stresses. The accumulation pattern of cellulose content under different abiotic stresses was similar to the expression trend of AsCesA genes. Our results provide valuable insights into the role of cellulose synthase in A.sinensis in plant defense.
ABSTRACT
ObjectiveTo clone coumarate-3-hydroxylase gene （C3H） from Angelica sinensis， and analyze the correlation between its bioinformatics， expression patterns and content of ferulic acid， and to explore the functions of ASC3H. MethodReal-time polymerase chain reaction （Real-time PCR） was used to clone the full-length cDNA of ASC3H based on the transcriptome dataset of A. sinensis， and the bioinformatics analysis of the gene sequence was carried out. Real-time PCR and high performance liquid chromatography （HPLC） were used to determine relative expression of ASC3H and content of ferulic acid in different root tissues of A. sinensis （periderm， cortex and stele）. ResultThe open reading frame （ORF） of ASC3H （GenBank accession number： MN2550298） was 1 530 bp， encoding 509 amino acids， with a theoretical molecular weight of 57.86 kDa and an isoelectric point of 8.36. It was a hydrophilic protein that was located in the chloroplast with multiple phosphorylation sites and a transmembrane region， and contained a conserved domain CGYDWPKGYGPIINVW_P450 （383-399 aa） in cytochrome P450. Multiple amino acid sequence alignment analysis showed that ASC3H had high similarity with C3H from other plants， especially Ammi majus in Umbelliferae. The Real-time PCR revealed that ASC3H had different expressions in periderm， cortex and stele tissues of A. sinensis roots. It was found from HPLC that the cortex tissues had the highest content of ferulic acid， and the stele tissues had the lowest. ConclusionASC3H was successfully cloned from A. sinensis， and its sequence characteristics were understood more clearly， suggesting that ASC3H might be involved in the ferulic acid biosynthesis pathway of A. sinensis. This paper provided a basis for further studying the functions of the gene and exploring the biosynthesis and regulation mechanism of ferulic acid in A. sinensis， while laying the foundation for the genetic improvement of A. sinensis．
ABSTRACT
Superoxide dismutase (SOD) is a key enzyme that scavenge superoxide anion free radical (O2·-) in vivo, and plays an important role in plant growth and development and stress. In this study, according to the genome and transcriptome data of Salvia miltiorrhizae, 9 SOD genes were identified and the expression patterns of SOD family genes were further analyzed, including 5 Cu/Zn-SOD, 2 Fe-SOD and 2 Mn-SOD. On the basis of proteomic analysis, combined with transcriptome data, one full-length cDNA of Mn-SOD gene, namely SmMSD2 was cloned from Salvia miltiorrhizae. The results of amino acid sequence alignment and phylogenetic analysis showed that SmMSD2 protein belongs to the manganese superoxide dismutase (Mn-SOD) subfamily, and SmMSD2 protein shares high sequence identity with the Mn-SOD proteins of various plants that all contain a C-terminal conserved metal-binding domain "DVWEHAYY". The prokaryotic expression vector pMAL-c2X-SmMSD2 was constructed and transformed into E. coli BL21 expressing strain, and the target recombinant protein was successfully induced and its enzymatic properties were analyzed. Spatiotemporal expression analysis showed that SmMSD2 gene was expressed in all tissues, indicating that SmMSD2 gene was constitutively expressed at a stable level. Real-time quantitative PCR indicated that drought (15% PEG6000), abscisic acid (ABA) and indole-3-acetic acid (IAA) could induce the expression of SmMSD2 gene, suggesting that SmMSD2 may be involved in the response of Salvia miltiorrhizae to abiotic stress such as drought, as well as the signaling pathways of phytohormone ABA and IAA. These results lay the foundation for further elucidating the involvement of superoxide dismutase in the stress response and accumulation of active components of Salvia miltiorrhiza.
ABSTRACT
ObjectiveTo identify the functions of the AP2/ERF family members in Pinellia ternata and promote the genetic improvement of P. ternata varieties. MethodWe identified and conducted a systematic bioinformatics analysis of the AP2/ERF family member genes in P. ternata based on the three generations of transcriptome data. Real-time polymerase Chain reaction （Real-time） PCR was employed to determine the expression pattern of AP2/ERF genes in different tissues and under different stress conditions. ResultA total of eight full-length AP2/ERF family members were identified from the transcriptome data， which were classified into three sub-gene families： AP2， ERF， and DREB. The deduced AP2/ERF proteins in P. ternata had the length of 251-512 aa， the theoretical pI of 5.29-11.72， the instability index of 45.90-82.41， subcellular localization in the nucleus， and conserved domains and motifs. AP2/ERF genes were expressed in different tissues of P. ternata， with high expression levels in the leaf. The stress response experiments showed that PtERF1 mainly responded to NaCl stress. The expression of PtERF2 and PtERF4 was significantly up-regulated under low temperature and polyethylene glycol （PEG）-simulated stress. PtERF3 responded to both low temperature and NaCl stress. The expression of PtERF5 was induced by high temperature， low temperature， NaCl and PEG stress. The expression of PtERF7 was up-regulated under high temperature， while that of PtERF8 under low temperature. ConclusionThe AP2/ERF genes in P. ternata can respond to stress and have the potential functions of regulating photosynthesis and improving root stress resistance.
ABSTRACT
Glutamate receptor-like (GLR) is an important class of Ca2+ channel proteins, playing important roles in plant growth and development as well as in response to biotic and abiotic stresses. In this paper, we performed genome-wide identification of banana GLR gene family based on banana genomic data. Moreover, we analyzed the basic physicochemical properties, gene structure, conserved motifs, promoter cis-acting elements, evolutionary relationships, and used real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) to verify the expression patterns of some GLR family members under low temperature of 4 ℃ and different hormone treatments. The results showed that there were 19 MaGLR family members in Musa acuminata, 16 MbGLR family members in Musa balbisiana and 14 MiGLR family members in Musa itinerans. Most of the members were stable proteins and had signal peptides, all of them had 3-6 transmembrane structures. Prediction of subcellular localization indicated that all of them were localized on the plasma membrane and irregularly distributed on the chromosome. Phylogenetic analysis revealed that banana GLRs could be divided into 3 subclades. The results of promoter cis-acting elements and transcription factor binding site prediction showed that there were multiple hormone- and stress-related response elements and 18 TFBS in banana GLR. RT-qPCR analysis showed that MaGLR1.1 and MaGLR3.5 responded positively to low temperature stress and were significantly expressed in abscisic acid/methyl jasmonate treatments. In conclusion, the results of this study suggest that GLR, a highly conserved family of ion channels, may play an important role in the growth and development process and stress resistance of banana.
Subject(s)
Musa/metabolism , Phylogeny , Abscisic Acid/metabolism , Temperature , Stress, Physiological/genetics , Hormones/metabolism , Gene Expression Regulation, Plant , Plant Proteins/metabolism , Gene Expression ProfilingABSTRACT
OBJECTIVE@#Flavonoids are the bioactive compounds in safflower (Carthamus tinctorius), in which chalcone synthase (CHS) is the first limiting enzyme. However, it is unclear that which chalcone synthase genes (CHSs) are participated in flavonoids biosynthesis in C. tinctorius. In this study, the CHSs in the molecular characterization and enzyme activities were investigated.@*METHODS@#Putative chalcone biosynthase genes were screened by the full-length transcriptome sequences data in C. tinctorius. Chalcone biosynthase genes in C. tinctorius (CtCHSs) were cloned from cDNA of flowers of C. tinctorius. The cloned gene sequences were analyzed by bioinformatics, and their expression patterns were analyzed by real-time PCR (RT-PCR). The protein of CtCHS in the development of flowers was detected by polyclonal antibody Western blot. A recombinant vector of CtCHS was constructed. The CtCHS recombinant protein was induced and purified to detect the enzyme reaction (catalyzing the reaction of p-coumaryl-CoA and malonyl-CoA to produce naringin chalcone). The reaction product was detected by HPLC and LC-MS.@*RESULTS@#Two full-length CtCHS genes were successfully cloned from the flowers of safflower (CtCHS1 and CtCHS3), with gene lengths of 1525 bp and 1358 bp, respectively. RT-PCR analysis showed that both genes were highly expressed in the flowers, but the expression of CtCHS1 was higher than that of CtCHS3 at each developmental stage of the flowers. WB analysis showed that only CtCHS1 protein could be detected at each developmental stage of the flowers. HPLC and LC-MS analyses showed that CtCHS1 could catalyze the conversion of p-coumaryl-CoA and malonyl-CoA substrates to naringin chalcone.@*CONCLUSION@#CtCHS1 is involved in the biosynthesis of naringin chalcone in safflower.
ABSTRACT
Chalcone isomerase is a key rate-limiting enzyme in the biosynthesis of flavonoids in higher plants, which determines the production of flavonoids in plants. In this study, RNA was extracted from different parts of Isatis indigotica and reverse-transcribed into cDNA. Specific primers with enzyme restriction sites were designed, and a chalcone isomerase gene was cloned from I. indigotica, named IiCHI. IiCHI was 756 bp in length, containing a complete open reading frame and encoding 251 amino acids. Homology analysis showed that IiCHI was closely related to CHI protein of Arabidopsis thaliana and had typical active sites of chalcone isomerase. Phylogenetic tree analysis showed that IiCHI was classified into type Ⅰ CHI clade. Recombinant prokaryotic expression vector pET28a-IiCHI was constructed and purified to obtain IiCHI recombinant protein. In vitro enzymatic analysis showed that the IiCHI protein could convert naringenin chalcone into naringenin, but could not catalyze the production of liquiritigenin by isoliquiritigenin. The results of real-time quantitative polymerase chain reaction(qPCR) showed that the expression level of IiCHI in the aboveground parts was higher than that in the underground parts and the expression level was the highest in the flowers of the aboveground parts, followed by leaves and stems, and no expression was observed in the roots and rhizomes of the underground parts. This study has confirmed the function of chalcone isomerase in I. indigotica and provided references for the biosynthesis of flavonoid components.
Subject(s)
Isatis/genetics , Plant Proteins/metabolism , Phylogeny , Arabidopsis/genetics , Flavonoids , Cloning, MolecularABSTRACT
GI (GIGANTEA) is one of the output key genes for circadian clock in the plant. The JrGI gene was cloned and its expression in different tissues was analyzed to facilitate the functional research of JrGI. RT-PCR (reverse transcription-polymerase chain reaction) was used to clone JrGI gene in present study. This gene was then analyzed by bioinformatics, subcellular localization and gene expression. The coding sequence (CDS) full length of JrGI gene was 3 516 bp, encoding 1 171 amino acids with a molecular mass of 128.60 kDa and a theoretical isoelectric point of 6.13. It was a hydrophilic protein. Phylogenetic analysis showed that JrGI of 'Xinxin 2' was highly homologous to GI of Populus euphratica. The result of subcellular localization showed that JrGI protein was located in nucleus. The JrGI, JrCO and JrFT genes in female flower buds undifferentiated and early differentiated of 'Xinxin 2' were analyzed by RT-qPCR (real-time quantitative PCR). The results showed that the expression of JrGI, JrCO and JrFT genes were the highest on morphological differentiation, implying the temporal and special regulation of JrGI in the differential process of female flower buds of'Xinxin 2'. In addition, RT-qPCR analysis showed that JrGI gene was expressed in all tissues examined, whereas the expression level in leaves was the highest. It is suggested that JrGI gene plays a key role in the development of walnut leaves.
Subject(s)
Juglans/genetics , Phylogeny , Plant Leaves , Cloning, Molecular , Gene Expression Regulation, Plant , Plant Proteins/metabolismABSTRACT
Dihydroflavonol 4-reductase (DFR) plays an essential role in the biosynthesis of anthocyanin and regulation of plant flower color. Based on the transcriptome data of Cistanche tubulosa (Schenk) Wight, a full-length cDNA sequence of CtDFR gene was cloned by reverse transcription-polymerase chain reaction (RT-PCR). CtDFR contains an open reading frame (ORF) of 1 263 bp which encodes 420 amino acids with a predicted molecular weight of 47.5 kDa. The sequence analysis showed that CtDFR contains a nicotinamide adenine dinucleotide phosphate (NADPH) binding domain and a specific substrate binding domain. The expression analysis indicated that CtDFR was highly expressed in red and purple flowers, and the relative expression levels were 4.04 and 19.37 times higher than those of white flowers, respectively. The recombinant CtDFR protein was expressed in E.coli BL21 (DE3) using vector pET-28a-CtDFR and was purified. In vitro enzyme activity analysis, CtDFR could reduce three types of dihydroflavonols including dihydrokaempferol, dihydroquercetin, and dihydromyricetin to leucopelargonidin, leucocyanidin and leucodelphinidin. Subcellular localization analysis showed that CtDFR was mainly localized in the cytoplasm. These results demonstrate that CtDFR plays an important role in regulation of flower color in C. tubulosa and make a valuable contribution for the further investigation on the regulation mechanism of C. tubulosa (Schenk) Wight flower color.
ABSTRACT
italic>Artemisia argyi (A. argyi) is a Chinese herbal medicine in China. The main active components are volatile oils, flavonoids, and other compounds, which have various pharmacological activities. Methoxylated flavonoids are the main active ingredients in A. argyi. Flavonoid O-methyltransferase (FOMT) is a key enzyme in the O-methylation of flavonoids. In order to further understand the function and characteristics of FOMT proteins, this paper carried out the whole genome mining and identification of FOMT genes in A. argyi and performed phylogenetic, chromosomal localization, gene sequence characterization, subcellular localization prediction, protein structure, gene structure analysis, and expression pattern analysis. The results showed that a total of 83 FOMT genes were identified in the genome of A. argyi. The phylogenetic tree shows that FOMT genes are divided into two subgroups, CCoAOMT (caffeoyl CoA O-methyltransferase) subfamily (32 genes) and COMT (caffeic acid O-methyltransferase) subfamily (51 genes). Gene sequence analysis showed that the number of amino acids encoded by FOMT was 70-734 aa, the molecular weight was 25 296.55-34 241.3 Da, and the isoelectric point was 4.51-9.99. Compared with 32 members of the CCoAOMT subfamily, nearly 1/3 of the 51 members of the COMT subfamily were hydrophobic proteins and 2/3 were hydrophilic proteins. Subcellular localization prediction showed that more than 80% of CCoAOMT subfamily members were located in the cytoplasm, and 96% of COMT subfamily members were located in the chloroplast. COMT subfamily members have more motifs than CCoAOMT subfamily members. The N-terminal motifs of COMT subfamily proteins are relatively variable, while the C-terminal motifs are relatively conserved. Expression pattern analysis showed that CCoAOMT subfamily members were mainly expressed in roots, while COMT members were mainly expressed in leaves. Some FOMTs showed the tissue expression specificity by real-time quantitative PCR analysis, especially in leaves. In this study, we identified and analyzed the FOMT gene family in A. argyi, and provided a theoretical basis for further research on the function of FOMTs and the biosynthesis of methylated flavonoids in A. argyi.
ABSTRACT
Glycogen synthase kinase 3/SHAGGY-like kinase (GSK3) proteins play important roles in regulating plant growth, development, and stress response. In order to reveal the characteristics of GSK family members in the medicinal plant Senna tora L., in this study, we conducted the identification and expression analyses of GSKs in S. tora based on its whole genome data, combined with bioinformatics and gene expression research methods. The results showed that a total of nine S. tora GSK genes were identified, all of which contained the GSK characteristic kinase domains. All members were distributed on six chromosomes, the encoding amino acid length ranged from 465 to 943 aa, the protein molecular weight was from 33.57 to 88.83 kDa, and the average isoelectric point was 8.2. The StoSKs were divided into four evolutionary branches, and the StoSKs in the same evolutionary branch shared the same exon/intron structure and conserved motifs. The expansion of the StoSKs gene family was mainly due to segment duplication events, and there were 17, 11, 8 and 7 pairs of collinear genes with Glycine max, Medicago truncatula, Arabidopsis thaliana and Oryza sativa, respectively. The promoter regions of StoSKs mostly contained responses elements related to stress stimulation, growth and development, and hormone induction. Transcriptome data analysis showed that StoSKs were expressed in different tissues, with the highest expression level in roots. Quantitative real-time PCR (qRT-PCR) analysis indicated that StoSKs in different evolutionary branches displayed a synergistic expression pattern response to light, and most of StoSKs could rapidly respond to NaCl stress with significantly up-regulated expression. All the results provide a basis for further analysis of the biological functions of the GSKs gene family in S. tora.
ABSTRACT
Cytochrome P450 (CYP450) is a kind of superfamily oxidase containing heme, which is distributed in various aerobic organisms. They are widely involved in the biosynthesis of terpenoids, alkaloids, flavonoids, fatty acids, etc. In this study, the full-length cDNA sequence of a P450 was cloned by reverse transcription-PCR (RT-PCR) and rapid amplification of cDNA ends (RACE) technology, with the specific primers that designed according to the sequence of a transcript annotated as P450 from the Aquilaria sinensis (Lour.) Gilg transcriptome database. The tissue expression and subcellular localization were also studied. The full-length cDNA of the cloned P450 gene is 1 920 bp, with 88 bp 5′-untranslated region (UTR), 344 bp 3′-UTR and a 21 bp polyA tail, and 1 488 bp open reading frame (ORF), encoding 495 amino acids. Sequence alignment revealed that the protein belonged to CYP71D family of cytochrome P450 family, and named AsCYP71D1. Tissue expression analysis indicated that AsCYP71D1 was mainly expressed in stem. Further subcellular localization of onion epidermis showed that AsCYP71D1 was expressed in cytoplasm, nucleus and cell membrane. This study will provide a foundation for further research on its function in agarwood sesquiterpene biosynthesis.
ABSTRACT
2-Oxoglutarate/Fe(II)-dependent dioxygenases (2-ODD) play an important role in plant primary and secondary metabolism. Based on the high-throughput sequencing platform Illumina NovaSeq 6000, the transcriptome of Salvia apiana Jepson was sequenced, and the obtained reads were de novo assembled. A total of 38 534 unigenes were obtained from the transcriptome. The assembled unigenes were annotated and 29 982 unigenes were given functional annotations. The 2-ODD genes were identified from the assembled S. apiana transcriptome database by bioinformatics methods, and the genes were analyzed, including the homology of the sequences, physicochemical characteristics, signal peptides, transmembrane domains, subcellular localization, secondary structure and tertiary structure, etc. The evolutionary relationships and the expression patterns of the identified 2-ODD genes were also analyzed. 39 full-length 2-ODD genes were identified from the transcriptome of S. apiana. The average length of these 2-ODD encoding proteins was 320 amino acids, the molecular weight was about 36.00 kDa, and most of them were hydrophilic proteins. Phylogenetic analysis divided these 2-ODD genes into several subfamilies. Gene expression analysis indicated that the 2-ODD genes were expressed in different parts of S. apiana, and the expression level of most genes was much higher in roots than that in leaves. This study can lay a foundation for further study of 2-ODD genes in S. apiana.
ABSTRACT
4-(Cytidine 5′-diphospho)-2-C-methyl-D-erythritol kinase (CMK) was one of the key enzymes in the methylerythritol-4-phosphate (MEP) pathway to generate terpenoids. In this study, based on the transcriptome data of Atractylodes lancea, the sequence of the CMK gene was cloned, named AlCMK (GenBank accession number OM283293). The results showed that AlCMK contains a 1 230 bp open reading frame (ORF) encoding 409 amino acids. The deduced protein had a theoretical molecular weight of 44 752.53 and an isoelectric point of 6.67. Transmembrane structure analysis showed that there was no transmembrane structure, and the secondary structure of AlCMK was predicted to be mainly composed of random coil. Homologous alignment revealed that AlCMK shared high sequence identity with the CMK proteins of Tanacetum cinerariifolium, Osmanthus fragrans, Eucommia ulmoides, Lonicera japonica and Salvia miltiorrhiza. Phylogenetic analysis indicated that AlCMK protein had the higher homology with CMK protein of Compositae. The pET-32a-AlCMK prokaryotic expression vector was constructed and a fusion protein with molecular mass of about 65 kDa was expressed in the E. coli BL21 (DE3). The qRT-PCR was used to analyze the expression pattern of AlCMK gene in different tissues and after MeJA treatment. Meanwhile, the enzyme activity was determined by ELISA kit. The results showed that AlCMK gene was tissue-expressed in different origins and its expression was induced by MeJA, and the results of the enzyme activity assay showed that the AlCMK enzyme activity in different regions was higher in the leaves. The subcellular localization showed that AlCMK was located in the chloroplast. This study provides a reference for further elucidating the biological function of AlCMK gene in terpenoid synthesis pathway in Atractylodes lancea.
ABSTRACT
WRKY, a class of conserved transcription factors in plants, plays important roles in plant growth, development and secondary metabolism. In the present study, 65 WRKY members were identified from de novo transcriptome sequencing data of three different tissues (root, stems and leaves) of Baphicacanthus cusia. BcWRKY proteins contained from 221 to 706 amino acids and the isoelectric point is from 4.68 to 9.68. Molecular weights range from 25 711.8 to 75 475 Da. The main secondary structures of BcWRKYs protein are random coil. A subcellular localization prediction indicated that the putative BcWRKY proteins were enriched in the nuclear region. Phylogenetic analysis showed that BcWRKYs could be categorized into three groups and five subgroups (Group IIa, Group IIb, Group IIc, Group IId and Group IIe) in Group II. Structural analysis found that all BcWRKY proteins contained a highly conserved motif WRKYGQK. Finally, the transcriptional profiles of ten BcWRKY genes highly expressed in root, stem and leaf tissues under abscisic acid (ABA), methyl jasmonate (MeJA), or salicylic acid (SA) treatment were systematically investigated using qRT-PCR analysis. Results showed that a total of ten BcWRKY genes were differentially expressed in response to ABA, MeJA, and SA treatment. This work would be provided a basis for further elucidating the molecular mechanism of WRKY transcription factors in the biosynthesis of indole alkaloids in B. cusia.
ABSTRACT
ObjectiveThe type 2C protein phosphatases （PP2C） are involved in numerous plant signal transduction pathways. They mainly participate in plant stress response and regulate second metabolites biosynthesis via negatively regulating MAPK signaling pathway. Herein，we were to identify and analyze PP2C （CsPP2C） gene family from hemp genome，in hope of providing comprehensive insights for studying CsPP2C function during the development of hemp. MethodMolecular Evolutionary Genetics Analysis （MAGA）-X was used to construct phylogenetic tree. Expert Protein Analysis System （ExPASy），WoLF PSORT，Multiple EM for Motif Elicitation （MEME），Batch Conserved Domain Search （Batch-CD-Search），PlantCare，and TBtools were used，respectively，to predict CsPP2C physicochemical properties，subcellular localization，conserved motifs，protein structure，cis-element in promoter and collinearity with Arabidopsis PP2C. Cannabis sativa transcriptome and Real-time polymerase chain reaction（Real-time PCR） were used to analyze and verify gene expressions，respectively. ResultFifty-two CsPP2C with conserved domains were identified from the entire genome of hemp，encoding proteins ranging from 244 to 1 089 aa in length and with molecular weights ranging from 26.76 to 122.53 kDa. Those genes were mainly distributed in the nucleus，cytoplasm and chloroplast. The 47 CsPP2C were divided into 10 subfamilies，and the remaining 5 were not clustered. Seven pairs of homologous genes between hemp and Arabidopsis thaliana were identified according to collinear analysis. The light-responsive elements and abscisic acid elements are most abundant in the prediction. The gene expression heat map showed varied expression pattern of CsPP2C in different tissues. Real-time PCR results of three CsPP2C were consistent with transcriptome data. Moreover，alternative splicing analysis showed that some CsPP2C had alternative-splicing genes during evolution. ConclusionWe predicted and analyzed CsPP2C gene family in genomic scale and showed that CsPP2C are involved in many biological processes，whereby provides foundation for CsPP2C functional study.
ABSTRACT
ObjectiveThe type 2C protein phosphatases （PP2C） are involved in numerous plant signal transduction pathways. They mainly participate in plant stress response and regulate second metabolites biosynthesis via negatively regulating MAPK signaling pathway. Herein，we were to identify and analyze PP2C （CsPP2C） gene family from hemp genome，in hope of providing comprehensive insights for studying CsPP2C function during the development of hemp. MethodMolecular Evolutionary Genetics Analysis （MAGA）-X was used to construct phylogenetic tree. Expert Protein Analysis System （ExPASy），WoLF PSORT，Multiple EM for Motif Elicitation （MEME），Batch Conserved Domain Search （Batch-CD-Search），PlantCare，and TBtools were used，respectively，to predict CsPP2C physicochemical properties，subcellular localization，conserved motifs，protein structure，cis-element in promoter and collinearity with Arabidopsis PP2C. Cannabis sativa transcriptome and Real-time polymerase chain reaction（Real-time PCR） were used to analyze and verify gene expressions，respectively. ResultFifty-two CsPP2C with conserved domains were identified from the entire genome of hemp，encoding proteins ranging from 244 to 1 089 aa in length and with molecular weights ranging from 26.76 to 122.53 kDa. Those genes were mainly distributed in the nucleus，cytoplasm and chloroplast. The 47 CsPP2C were divided into 10 subfamilies，and the remaining 5 were not clustered. Seven pairs of homologous genes between hemp and Arabidopsis thaliana were identified according to collinear analysis. The light-responsive elements and abscisic acid elements are most abundant in the prediction. The gene expression heat map showed varied expression pattern of CsPP2C in different tissues. Real-time PCR results of three CsPP2C were consistent with transcriptome data. Moreover，alternative splicing analysis showed that some CsPP2C had alternative-splicing genes during evolution. ConclusionWe predicted and analyzed CsPP2C gene family in genomic scale and showed that CsPP2C are involved in many biological processes，whereby provides foundation for CsPP2C functional study.
ABSTRACT
MYB transcription factors, one of the largest transcription factor families in plants, play an important role in signal transduction, plant growth and plant resistance. In this study a full-length cDNA of the PnMYB1R1 gene was cloned from Panax notoginseng. Sequence analysis, prokaryotic expression and purification, subcellular location, transcriptional activity analysis, tissue-specific analysis and expression analysis under different abiotic stresses was performed. The open reading frame (ORF) of PnMYB1R gene was 738 bp, encoding a protein of 245 amino acids with a predicted molecular mass (MW) of 27.0 kD. The sequence analysis and polygenetic analysis indicated that the PnMYB1R1 protein contains a conserved R3 domain, belonging to TRF-like protein in 1R-MYB-type transcription factors. The recombinant PnMYB1R1 protein was expressed in Escherichia coli BL21(DE3) cells using the prokaryotic expression vector pET28a-PnMYB1R1 and was purified. Subcellular localization analysis showed that PnMYB1R1 was localized in the nucleus. Transcriptional activity analysis indicated that the PnMYB1R1 transcription factor has transcriptional activation activity. Expression analysis indicated that PnMYB1R1 was primarily expressed in roots, followed by stems and leaves, and then rootlets. The expression level of PnMYB1R1 in root, stems, leaves and rootlets was influenced by salt, low temperature and drought treatment, while the abundance of PnMYB1R1 was significantly induced by salt stress in these tissues. These results provide valuable insights into the role of 1R-MYB transcription factors in plant defense.
ABSTRACT
Licorice, one of the most commonly used medicinal materials in China, grows mainly in arid and semi-arid regions and has important economic and ecological values. Basic leucine zipper (bZIP) transcription factors in plants play an important role in regulating biological or abiotic stress responses, growth, and secondary metabolite synthesis. bZIP transcription factors in the published whole genome database of Glycyrrhiza uralensis were identified using bZIP sequences found in Arabidopsis thaliana genome as reference, and ABA-dependent bZIP genes were identified by using Illumina high-throughput sequencing. The physical and chemical properties, structure of the encoded proteins, and the gene expression patterns with exogenous ABA stress were analyzed. A total of 69 bZIP transcription factor genes were identified in G. uralensis, named Gubzip1-69, and they were divided into 10 subfamilies (A-I and S) according to their similarity to bZIPs of A. thaliana. By calculating the relative expression levels of the 69 GubZIPs genes under different concentrations of exogenous ABA stress, genes that may be involved in the regulation of ABA signaling pathways were identified, namely GubZIP1, GubZIP5, GubZIP8, GubZIP30, GubZIP33 and GubZIP56. The results of expression pattern analysis of these GubZIPs genes under exogenous ABA stress showed that the expression pattern of GubZIPs genes changed significantly with 50 mg·L-1 ABA. The relative expression levels of these genes decreased 3 h after treatment, and gradually increased 6 h after treatment. Except for GubZIP8, the relative expression levels of these genes were significantly increased after 12 h. Further research on the function of bZIP transcription factors of G. uralensis and elucidating their regulatory mechanisms should be of interest and will provide a scientific basis for cultivating high-quality cultivars of G. uralensis through molecular breeding methods.
ABSTRACT
A large number of putative risk genes for autism spectrum disorder (ASD) have been reported. The functions of most of these susceptibility genes in developing brains remain unknown, and causal relationships between their variation and autism traits have not been established. The aim of this study was to predict putative risk genes at the whole-genome level based on the analysis of gene co-expression with a group of high-confidence ASD risk genes (hcASDs). The results showed that three gene features - gene size, mRNA abundance, and guanine-cytosine content - affect the genome-wide co-expression profiles of hcASDs. To circumvent the interference of these features in gene co-expression analysis, we developed a method to determine whether a gene is significantly co-expressed with hcASDs by statistically comparing the co-expression profile of this gene with hcASDs to that of this gene with permuted gene sets of feature-matched genes. This method is referred to as "matched-gene co-expression analysis" (MGCA). With MGCA, we demonstrated the convergence in developmental expression profiles of hcASDs and improved the efficacy of risk gene prediction. The results of analysis of two recently-reported ASD candidate genes, CDH11 and CDH9, suggested the involvement of CDH11, but not CDH9, in ASD. Consistent with this prediction, behavioral studies showed that Cdh11-null mice, but not Cdh9-null mice, have multiple autism-like behavioral alterations. This study highlights the power of MGCA in revealing ASD-associated genes and the potential role of CDH11 in ASD.